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Abstract 

We define and investigate the properties of the jaggedness of path integral trajecto- 
ries. The new quantity is shown to be scale invariant and to satisfy a self-averaging 
property. Jaggedness allows for a classification of path integral trajectories accord- 
ing to their relevance. We show that in the continuum limit the only paths that 
are not of measure zero are those with jaggedness 1/2, i.e. belonging to the same 
equivalence class as random walks. The set of relevant trajectories is thus narrowed 
down to a specific subset of non-differentiable paths. For numerical calculations, we 
show that jaggedness represents an important practical criterion for assessing the 
quality of trajectory generating algorithms. We illustrate the obtained results with 
Monte Carlo simulations of several different models. 
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1 Introduction 



The set of all the trajectories that one integrates over in a path integral is huge 
(cardinality K2). In order to have successful numerical simulations it is impor- 
tant to generate a finite number of trajectories that are as representative of the 
whole set as possible. For this reason it is necessary to broaden the amount of 
analytical information regarding the subset of trajectories that give dominant 
contributions to path integrals. This subset is much smaller than the set of all 
trajectories. Investigation of its properties is of great importance. Ultimately 
we would like to find ways to generate only these relevant trajectories. 
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Classification of trajectories according to their relevance is just as important 
for analytical calculations of continuum limit amplitudes. As a result of the 
continuum limit many trajectories are in fact of measure zero, i.e. give no con- 
tribution to path integrals. For example, it is well known that differentiable 
trajectories are of measure zero. This has been a severe impediment to the 
development of path integration theory as most of mathematics (and practi- 
cally all of our intuition) is firmly based on smooth, differentiable functions. 
In fact, the principle reason that we have made any analytical progress with 
path integrals lies in the fact that, although of measure zero, some differen- 
tiable functions (e.g. classical solutions) turn out to be important markers of 
"nearby" non-differentiable paths that do give important contributions. 

We are in a precarious position in which we know very little about the paths 
we need to work with. In addition, most of our knowledge is negative, i.e. 
tells us which trajectories are not important rather than which are. This is a 
substantial problem in the analytical approach. In numerical simulations, on 
the other hand, we work with discrete trajectories about whose classification 
we know even less. Therefore, the implementing of a systematic classifica- 
tion of the relevance of path integral trajectories (both in the discretized and 
continuum theories) is a crucial starting point in the quest for more efficient 
calculation schemes. 

One of the few positive statements concerning path integral trajectories is that 
they are stochastically self-similar [1,2]. An exhaustive review of various as- 
pects of the path integral method can be found in [3] . A direct consequence of 
self-similarity is that relevant path integral trajectories have [4] fractal (Haus- 
dorff) dimension du = 2. In this paper we introduce a new classification 
property - the jaggedness of a trajectory. An analytical and numerical analy- 
sis of the properties of jaggedness leads to a new classification of trajectories 
according to their relevance to path integrals. 



2 Self-similarity of trajectories 

The property of stochastic self-similarity of path integral trajectories has im- 
portant repercussions both on the dynamics as well as on the construction of 
efficient numerical algorithms for generating paths. In a previous set of pa- 
pers [5,6,7] we have used self-similarity to obtain analytical relations between 
discretizations of different coarseness for the case of a general theory. The 
newly developed analytical method systematically improves the convergence 
of path integrals of a generic iV-fold discretized theory through the explicit 
construction of a set of effective actions for p — 1, 2, 3, . . .. These effective 
actions are equivalent to the starting action (in the sense that they lead to 
the same continuum amplitudes), however, the path integrals calculated us- 
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ing them converge to the continuum limit ever faster. Discretized amplitudes 
calculated using the p level effective action tend to the continuum limit as 
1/N P . Using the general procedure we obtained explicit effective actions up 
to p — 9. Self-similarity played a crucial role in this procedure allowing us 
to derive, and asymptotically solve, an integral equation relating discretized 
theories viewed at different coarseness. It was shown [7] that this approach is 
in fact equivalent to the derivation of a generalization of Euler's summation 
formula to path integrals. 

Self-similarity has also been successfully utilized in constructing efficient nu- 
merical algorithms for generating paths for Monte Carlo simulations [8,9]. In 
particular the Levy construction discussed in these references generates self- 
similar paths through a simple iterative procedure. One begins with the fixed 
end points a at time t — and b at time t = T and samples a bisecting point 
at time t — T/2 from a Gaussian centered in the middle between a and b and 
with width a = \Jt /2. Having sampled the point at t = T/2 one now bisects 
the two new intervals [0, T/2] and [T/2, T] in the same way generating new 
bisection points at t — T/4 and t = 3T/4 with appropriately centered Gaus- 



sians of width a = JT/A. The procedure is continued recursively doubling 



the number of sampled points at each level and using the width a = v/e/2 



where e is the current time step. This method exactly samples free particle 
trajectories, but also works very well for interacting theories. The principle 
benefit of the method is that computational effort scales as O(N), where N 
is the coarseness of the discretization. 

It is important to strengthen the relation between these two types of ap- 
proaches. For a deeper understanding of the role of self-similarity it is nec- 
essary to classify trajectories according to their relative contribution to the 
path integral. As we have already mentioned, the first and most natural such 
classification of paths was with respect to their fractal dimension. The fractal 
dimension dn is calculated from the formula 



where (L) is the expectation value of L = X^q 1 \q i+1 — qi\ : the total trajectory 
length. It was shown by Kroger [4] that for theories with (Euclidean) action 
of the form 



the only trajectories that contribute to the path integral are those with fractal 
dimension dn — 2. The reason for this is that for short times of propagation 
the kinetic term dominates over the potential and each model looks like a 
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random walk (which has cLh = 2). This is illustrated in Fig. 1 on the example 
of an anharmonic oscillator with quartic coupling V(q) = |g 2 + -^gq 4 - Kroger 
has also shown that the addition of velocity dependent interactions changes 
the fractal dimension of relevant paths. Let us now introduce another quantity 
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Fig. 1. (Left) A typical graph of the expectation value of the trajectory length (L) 
as a function of discretization coarseness N. The data fits to (L) = 0.8y/N, i.e. 
to dfj = 2. (Right) A detailed picture of how the fractal dimension approaches the 
value for a random walk in the continuum limit fits to 2 — dn = 0.0A/y/N. Both plots 
correspond to propagation from a = to b = 1 in time T = 1 for an anharmonic 
oscillator with quartic coupling g = 1. The number of Monte Carlo samples used 
was N M C = 9.2 • 10 6 . 




Fig. 2. An example of a trajectory that goes from a (at t = 0) to b (at t = T) in 
N = 14 discrete time steps. 

(complementary to the fractal dimension) which will serve to further classify 
path integral trajectories according to their relevance. 



The iV-fold discretized path from q = a to qN = b is determined by the N — 1 
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intermediate positions qi, q 2 , . . . , qN-i- We define the jaggedness of a path as 
1 N ~ 2 1 

J= W —rJ2^-^m + i)) , (3) 
iV 1 i=0 z 

where 5i = — From the above definition we see that J in fact counts the 
number of peaks (both minima and maxima) divided by N — 1 (the maximal 
number peaks for an A-fold discretized trajectory). Therefore, for all values 
of N, we have J G [0, 1]. Fig. 2 illustrates a typical path with N = 14 and 
jaggedness J = -j§ (intermediate points that are peaks are depicted by black 
circles). 

An important property of the jaggedness is that it is scale invariant up to 
1/N terms. Namely, if we scale the intermediate points q^ — > \q,i the only 
thing that can change are the first and last term in the above sum (since 
the end points need to remain fixed). This brings about a change in J that 
is 0(1/N). Therefore, we see that in the continuum limit J is in fact scale 
invariant. J is also scale invariant for finite N when a = b = 0. This behavior 
is quite different from that of fractal dimension. From eq. (1) we see that the 
fractal dimension is not invariant under scaling with an A-dependent A. For 
this reason the classification of trajectories with respect to J is independent 
of a classification with respect to dn- 



3 Analytical and numerical analysis of jaggedness 

In the continuum limit all smooth and differentiable trajectories have J = 0. 
This is not only true of monotonic trajectories but also of those having a finite 
number of extrema. Paths with J = are of measure zero. A classification 
of paths with non- vanishing jaggedness is therefore a path integral motivated 
classification of non-differentiable paths. 

From eq. (3) we see that (up to 1/N terms) the jaggedness satisfies an aver- 
aging property, i.e. if we split a trajectory with 2N discrete time steps into 
two equal halves (with jaggedness J\ and J 2 ), the total jaggedness equals 

1 - sgn(£jv„i£jy) _ Ji + J 2 

2 2 J U 

As we can see, for finite N the averaging is broken by a 1/N term. In the 
continuum limit, however, the averaging of jaggedness is exact. Cutting up a 
trajectory into k equal pieces we get J = \ Yh=i Ji- Stochastic self-similarity 
now implies that all the Jj's are equal to each other, and in fact that they are 



J = -(Ji + J 2 ) 

2 V ; 2N- 1 
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equal to the jaggedness of the whole path. We denote this property of jagged- 
ness as self-averaging. The jaggedness of the whole trajectory (corresponding 
to a motion for time T) is thus equal to the jaggedness of even the smallest 
piece of that trajectory (corresponding to the propagation for a much shorter 
time At). 

It is well known that for short times of propagation the dynamics of any model 
is well approximated by that of a random walk, i.e. does not depend on the 
potential. Therefore, self-similarity of trajectories implies that the jaggedness 
of the trajectories of a general model with action of the form given in eq. (2) 
is equal to the jaggedness of trajectories for a random walk. 

The expectation value of the jaggedness for a random walk is quite easily 
calculated. For a random walk S{ and 6{+i are not correlated and so (J) = |. 
We therefore conclude that for all the models with action of the form of eq. 
(2) we have (J) = |. Similarly, for finite N we find that (J) deviates from the 
continuum as 0(1/ N). 
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Fig. 3. Relative contributions to the path integral of trajectories of different jagged- 
ness for different models. The histograms are practically indistinguishable for: an- 
harmonic oscillators with quartic coupling g = 1 (black) and g = 1000 (red) as well 
as particles in a modified Poschl- Teller potential V(q) = with param- 

eters a = 0.5, (3 = 1.5 (green) and a = 0.5, j3 = 2 (dashed black). The data is for 
N = 512, T = 1, a = 0, b = 1, Nmc = 9-2 • 10 6 and fits to a Gaussian centered at 
0.5 with width 0.022. The number of bins is 100. 

Fig. 3 gives the relative contributions to the path integral of trajectories of 
different jaggedness for the case of N = 512. As is shown in the figure, the 
distribution is practically indistinguishable for different models and widely 
varying parameters. Detailed numerical investigations show that for all models 
studied deviations from a Gaussian distribution (e.g. skew and kurtosis) go to 
zero as N — > oo. In agreement with our previous analytical argument based on 
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Fig. 4. Width of the distribution around (J) plotted as a function of N. The data 
fits to a j = 0.5 /y/N. The plot is for an anharmonic oscillator with quartic coupling 
g = 1, propagating for T = 1 from a = to b = 1 and iVwc = 9.2 • 10 6 . 

self-similarity the center of the distribution is (J) = | + 0(1/N). The width 
of the distribution oj vanishes as oj ~ 1/y/N as illustrated in Fig. 4. The 
conclusion is that for models with action of the form of eq. (2) the only paths 
that are not of measure zero are those with J = |. 

There are several important consequences of this that we need to mention. 
The analytical consequence of the above results is that they give us a much 
more detailed understanding of which trajectories are relevant and which are 
not. Not only are smooth, differentiable trajectories of measure zero (i.e. those 
with J = 0), but in fact most of the non-differentiable trajectories also do not 
contribute to path integrals. To understand path integrals better we need to 
focus on a much narrower class of non-differentiable trajectories - those with 
J = |, i.e. those belonging to the same jaggedness equivalence class as the 
random walk. 

The numerical consequence of the above results is that for finite N it is im- 
portant to have algorithms that generate trajectories that are distributed in 
a way that mimics, as much as possible, the physical distribution of relative 
contributions to the path integral such as the one shown in Fig. 3. The most 
important thing to focus on is the center of the distributions. 

Fig. 5 compares the physical average (J), which is algorithm independent, 
with J the average jaggedness of the trajectories generated by the algorithm 
(in this case the Levy construction). Note that for the Levy construction J is 
completely independent of the dynamics. The reason why this method repre- 
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Fig. 5. iV-dependances of (J) and of J. The fact that J tends to the same continuum 
value as the physical expectation value indicates that the Levy construction is a good 
algorithm to use for large N. Parameters are the same as in Fig. 4. 

sents a good algorithm (at least for large N) is that J — > ^ in the continuum 
limit. This figure also indicates that a useful measure of the quality Q of an 
algorithm for generating trajectories is roughly the inverse of \(J) — J\. We 
will use this measure of quality to compare the Levy construction with the 
diagonalization algorithm in which trajectories are generated by a Gaussian 
distribution function using a semi-classical expansion. The computing time of 
this algorithm scales as 0(N 2 ) since it is necessary to diagonalize the quadratic 
form in the exponential of the distribution function. 

In Fig. 6 we plot |(J) — J| for these two algorithms. The top curve is for 
the Levy construction, the bottom for the diagonalization method. The latter 
algorithm outperforms the former by an order of magnitude for all values of 
N. As a result of this the Monte Carlo errors calculated using the two methods 
differ by an order of magnitude. However, one should not forget that quality 
needs to be balanced by an assessment of cost (in computing time) . Although 
it gives a relatively lower quality, the computing time for the Levy method 
scales as N. The diagonalization method outperforms the Levy method for a 
given N, however its computation time scales as N 2 . Both algorithms are good 
(since \ (J) — J\ vanishes in the continuum limit). A simple cost-benefit analysis 
tells us that Levy construction is better for large N, while diagonalization is 
better for smaller N. 

Not all algorithms for generating paths are good, however. For example, if 
we modify the Levy construction using uniform distributions (of appropriate 
widths) instead of Gaussians we get an algorithm for which J does not tend 
to |. A signature of a bad algorithm is that there exists a maximum quality 
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Fig. 6. Comparison of (J) and J calculated using two different algorithms for the 
case of an anharmonic oscillator with quartic coupling. The top curve was calculated 
using the Levy bisection method. The data fits to the curve \ (J) —J\Levy = 0.051/iV. 
The lower curve was calculated using the diagonalization method. The data fits to 

\(J)-J\diag = 0.003/N. 

Qmax that can't be passed irrespective of computational cost. 

As we have mentioned, it is important to strengthen the connection between 
the dynamical (calculation of path integrals) and kinematical (generation of 
trajectories) aspects of self-similarity. The effective action approach [5,6,7] 
has brought about an immense speedup in path integral calculations. The 
speedup is a direct consequence of the fact that effective actions allow us 
to work with much smaller values of iV to obtain the same precision. Using 
standard algorithms (which have all been derived to be optimal for large N) 
we have obtained a speedup of many orders of magnitude. The important step 
that next needs to be taken is to develop a path generating algorithm that is 
tailored for small iV's, i.e. for which J is near to (J) for coarse discretizations. 
Note that the diagonalization algorithm is one such method, but that it is 
not computationally optimized. We are currently working on developing these 
kinds of algorithms and tying them in to the effective action approach. 

At the very end we wish to make further contact between jaggedness and 
random walks. In order to be near the random walk limit, the potential term 
in the discretized action eV, where e = T/N, must be smaller than the kinetic 
term 5 2 /2e. One should expect that when we are near the random walk regime 
all quantities depend on the ratio of these two. On the other hand, for a 
random walk we have 5 2 /e ~ 1, so that the ratio is in fact eV(q c ) where q c is a 
characteristic length. A rough value for the characteristic length follows from 
q c p ~ 1 (essentially Heisenberg's uncertainty relation in h = 1 units). We are 
using units in which m = 1, so p = 5/e ~ 1/ sfe. The last step follows from the 
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basic random walk relation 5 2 /e ~ 1. Finally we find that everything should 
be expressed in terms of the ratio eV(yfe). For example, for the anharmonic 
oscillator with quartic coupling this ratio is e 3 g. We have shown that (J) 
differs from | by a term proportional to 1/N. From this it follows that for 
the oscillator with quartic anharmonicity \(J) — || should be proportional to 
g 1 / 3 /N. A similar back of the envelope calculation for a particle moving in a 
modified Poschl- Teller potential gives that |(J) — -A should be proportional to 

- Fig. 7 illustrates that these simple calculations do in fact hold. 




Fig. 7. (Left) Anharmonic oscillators with quartic coupling g = 1, 2,4, 8, . . . , 64 for 
N = 128,256 and 512. (Right) Modified Poschl-Teller potential for a = 0.5 and 
(3 = 1.5,2.5,... ,6.5 for N = 128,256 and 512. All the simulations were done for 
a = 0, b = 1, T = 1 and N MC = 9.2 • 10 6 . 

To conclude, we have identified and investigated the properties of a quantity 
that we call the jaggedness and that is useful for obtaining a more detailed 
classification of relevant path integral trajectories. For discrete calculation, i.e. 
numerical simulations, the properties of the jaggedness are useful for obtaining 
more efficient algorithms for generating representative paths. Furthermore, we 
have shown that jaggedness can be used as an important practical criterion 
of the quality of trajectory generating algorithms. In the continuum limit 
(analytical calculations) we found that only trajectories with jaggedness equal 
to \ contribute to the path integral. In this way we greatly narrow the set of 
trajectories that are not of measure zero to those belonging to the equivalence 
class of the random walk. Classification of paths with respect to jaggedness 
is thus a classification of relevant non-differentiable paths. In the continuum 
limit, jaggedness is shown to be scale invariant as well as self-averaging. 
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